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The overall intent of this work was to investigate the physics of high energy 
photons and particles (especially electrons and positrons) in the compact, high-energy- 
density environment of active galactic nuclei and quasars (”AGN”). These objects are 
observed to contain hard X-ray and 7-ray (with hu^m = ec *) photons and relativistic 
electrons and/or positrons (with E » m e c ^). The high energy density, and consequent 
high optical depth, of the nuclear region means these particles are very likely to interact 
with each other as well as with their surroundings before emerging from the nucleus 
and being observed. 

The fundamental goal of the work was to determine the effect of this compact, 
optically thick ’’atmosphere” on the emergent particle and photon spectra. This was 
pursued by developing and applying a numerical code which follows the nonlinear 
spect”ral evolution of a pair/photon plasma, due to two-body scattering and inter- 
action processes, in an unmagnetized system. The code was applied both to static 
plasmas and to relativistic, expanding winds. The first application, while not physi- 
cally correct for long times (when other processes, not included in the code, must come 
into play), explored the fundamental effects of the microphysics on the spectral evolu- 
tion. The second application is probably the relevant one for active galactic nuclei: the 
expansion of an internally relativistic plasma will at some point shut off ("freeze out”) 
the reactions and determine what an external observer sees. 

Several related tasks were performed during the grant period: 

a) developing a formalism to account for the effect of mass creation and annihilation 
on a relativistic wind; 

b) modelling coherent radio emission from relativistic electron beams in active galac- 
tic nuclei; 

c) writing and presenting an invited review on the role of high energy pairs and 
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photons in active galactic nuclei, at the ’’Supermassive Black Hole” meeting held 
to George Mason University, October, 1986. 

This work is summarized in this Final Report. 


Code Development and Cross Sections 

The overall picture assumed in the calculation is that of an optically think, 
isotropic pair/photon fireball, with no escape of pairs of photons, nor any radiative 
transfer. This might, for instance, represent the optically thick atmosphere of the hot 
inner regions of a thin accretion disk, or the nonlinear development of an electrody- 
namic pair/photon cascade in a nuclear Poynting flux beam. 

The evolution of the pair or photon spectrum - n a (E ) - due to two-body processes 
is given by 

= E / ME’)n s (E")v^^(E\ E")g^(E\ E"-, E)dE'dE" (1) 

0,6 

where apjj is the total cross section for species 0 and 6 at energies E' and E"\ vpj is the 
relative velocity of the two quanta; and gp jj is the distribution of output energy E from 
this reaction. The indices a , 0 and 8 take the values e (for pairs) or 7 (for photons). 
Equation (1) was solved numerically by discretizing energy space in logarithmic bins, 
so that the number in bin i is 


n a i = [ n a (E)dE (2) 

Jbin i 

and the reaction rate coefficient is 

= HE>i3j,6k°il,s(EvE l ')gi3, s {Ej,E k -,E)dE (3) 

ca T Jbin i 

The distribution within the bin, f(E), was assumed to be f{E) oc 1 /E. If equation (1) 
is normalized to some fiducial density, n 0 , and to a characteristic collision time, using 
t = cajTiot, the discretized equation becomes 


dr 


0,6 j,k 


( 4 ) 


- 3 - 



The equations represented by (4) were solved simultaneously, numerically, to get the 
spectral evolution, n m -(r), for a = e, 7 and for each i. 

The differential cross sections, ^ necessary for this calculation were evaluated 
for four two-body reactions: pair annihilation, two-photon pair creation, Compton scat- 
tering, and electron-electron bremsstrahlung; these are the most important reactions 
at moderate energies (~ MeV to GeV) in unmagnetized plasmas. Calculating ^ 
required integrating the full, differential cross sections - which are available in the liter- 
ature - over scattering angles (in the center of momentum or electron rest frames) and 
enegies consistent with bin i. This was done numerically by Monte Carlo methods. The 
energy space bins ran from 3xlQ~^m e c* to 30m e c^, in steps of Alog(£^/m e c^) = 0.5. 
Imposed ” rebinning” of the output of the Monte Carlo code ensured energy and number 
conservation to very good accuracy at each step. 

2. Evolution of Non-Expanding systems. 

This code was developed to model expanding fireballs, as in active nuclei - where 
the conditions axe such that the two-body processes listed above will be important, 
but where expansion of the fireball will freeze out the spectrum before other processes 
(such as three-body processes, free-free absorption, etc) have had time to occur. We 
began, however, by using the code to follow the evolution of a confined, constant- 
density pair/photon plasma under the four processes listed above. This allowed the 
effect of the reactions to be studied in isolation, separate from the effects of the fireball 
expansion. Clearly, this code cannot follow the evolution of a closed system toward 
true thermal equilibrium, as inverse and three-body processes are not included. This 
study will, however, provide a good description of the early phases of thermalization, 
when the two-body processes included are the most important ones, and thus will allow 
intelligent interpretation of the results of the expanding- wind calculations. 

The numerical runs started with an arbitrary initial pair and/or photon spec- 
trum - the ” creation” spectrum. This might, for instance, be envisioned to be the 
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pair/photon spectrum produced by a hot, inner accretion disk around a massive black 
hole; or particle pickup and acceleration in a Poynting flux beam. Four types of ” cre- 
ation” spectra were chosen - monoenergetic pair and photon spectra, and ’’white” pair 
and photon spectra (in which n(E ) oc E~^, giving a flat energy spectrum). About 
fifteen different ’’creation” spectra at various energies within these categories were fol- 
lowed to long times. 

It was found that each experiment evolved as follows. The early phases of the 
evolution consisted of a slow evolution away from the initial spectrum, toward a com- 
mon, ’’quasi-equilibrim” spectrum. This spectrum is maintained by non-local (in energy 
space) scattering and energy transfer processes; thus it reflects the two-body differential 
cross sections rather than the initial conditions. Both pairs and photons share this com- 
mon shape at most energies; the pair/photon ratio reflects the creation/ annihilation 
rate balance. This common spectral shape was reached initially at moderate (~ m e c 2 ) 
energies. On a longer timescale, higher and lower energies reached this common spec- 
trum. In particular, initial high-energy (>> m e c 2 ) quanta were eventually degraded by 
interactions with cooler quanta; low-energy quanta were created by bremsstrahlung and 
by scattering. Two unexpected and important trends appeared in these experiments. 

The first interesting result concerns the rate at which the system reaches this 
mid-energy ’’thermal” peak. This rate might, naively, be expected to be controlled by 
the dominant reaction likely in the creation spectrum - two-photon pair creation in 
a photon spectrum, and pair annihilation in a pair spectrum. Indeed, such "linear” 
esrtimates have been comonly used in the literature. However, the strongly non-linear 
nature of the system does not bear out this approach. Monoenergetic starts were 

e\ 

found to evolve much faster than the linear estimates would predict (by a factor ~ 10 , 
typically; the exact rate depended on the specific creation spectrum). White starts 
were found to evolve much more slowly than the linear estimate (again by a factor 
£ 10 2 ). 
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The second interesting result is the emergence of the slowly evolving, quasi- 
equilibrium spectrum. This spectrum appears on timescales that should just be reached 
in the most compact AGN, before the wind expansion shuts off the evolution. This 
quasi-equilibrium spectrum has a broad peak, centered initially at E ~ 0.1 m e c 2 , in 
both n e (E) and riy(E). This peak decays only slowly to higher energies, being re- 
duced by ~ 10 3 at ~ 30 E for instance. The spectrum at the highest initial energies 
slowly degrades from its initial state to an approximate exponential tail. The low- 
energy side of the broad peak has n(E) oc if 1 / 2 for both pairs and photons; for pairs 
this continues to the lowest energy in the calculation, while for photons the spectrum 
at the lowest energies is approximately fiat, reflecting the creation of soft photons 
by bremsstrahlung. The form n(E) oc E 1 ^ is the low-energy limit of the Maxwell- 
Boltzmann distribution. This quasi- equilibrium evolves slowly, at a rate determined 
by the rate of bremsstrahlung creation of photons by pairs at the mean energy, E. As 
this process continues, E decays and the photon creation rate slows down even further. 
While the average energy decays, both the pair and photon spectra keep their charac- 
teristic shape, but the ratio n e (E)/n- r (E) decays as fewer photons exist with energies 
above the pair creation threshold. 

None of the runs was followed much beyond this point, both because other pro- 
cesses not included in the code would become important, and because the long-term 
evolution is probably not important in AGN. 

3. Incorporating the Relativistic Wind 

In the AGN environment, a relativistic pair/photon plasma with E > GMm e /r 
(and with GMm e /r ;$ m e c 2 ) is unlikely to remain trapped; rather a relativistic expan- 
sion is likely. Such an expansion will modify the spectral evolution described above by 
’’freezing out” the reactions at some intermediate point. The emergent spectrum will 
probably not be the fully thermalized, quasi-equilibrium described above, but rather 
some intermediate shape reflecting the effective optical depth of the flow and, for ’’thin 
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flows”, the creation spectrum. 


In order to include these effects in the models, one must know the velocity field 
of the wind, v(r ) = cj 3 (r ) (assuming a steady flow), or the bulk Lorentz factor, 7 (r) = 
(1 — . When the internal energy of the fluid is well above the gravitational 

potential, the flow is determined by 


dj 


dp 


W d? = ~ T dr- 


(5) 


where w = + pc^ is the enthalpy of the fluid; p is the density and p is the pressure 

of the fluid; T = dlnp/dlu. p is the adiabatic exponent. If mass were conserved, this 
equation could be reduced to one in only 7 (r) by adding a continuity equation and an 
equation of state (usually assumed to be adiabatic, p/p^ = constant). However, pair 
creation and annihilation means the adiabatic exponent is a complicated function of 
local conditions, and cannot easily be incorporated into (o). 

Instead, dp /dr was evaluated explicitly at each step. This was done by noting 
that, in terms of the binned n;, 



"7" 


( 6 ) 


(noting that E{ = e 2 m e c^ is the totcil. not kinetic, energy of the ith bin). The pressure 
gradient in (5) can then be evaluated explicitly by chain rule, with a pseudo- continuity 
equation invoked to account for the rate of change of specific volume in the flow. This 
leads to 


dp 

dr 


1 ( 2 

'3 


7 


d 7 


7 


2 — 1 dr 


A(r) + R(r) 


( 7 ) 



! 


and 


1 


R(r) = -m e c 2 J2 


dn 


ei 


e i I + 

ei 


dn 




dr V ~ l ei ) ' dr 

and dn e {/dr and dn-^/dr are the rates of change due to the two-body reactions; c.f. 
equation (4), and equation (9), below. 

Combination of (7) with (5) allows explicit evaluation of 7 (r) within the calcula- 
tion. 


4. Spectral Evolution in an Expanding Wind 

The case of an expanding wind was modelled by following the evolution of the 
flux at radius r, 

/ = 4tt r 2 nj3c/F 0 ( 8 ) 

(normalized to some constant flux, F 0 ). The number density in ( 8 ) is that in the 
observers frame; in the co-moving (fluid) frame, n? = n/ 7 . The binned rate equation, 
analogous to (4), becomes 


df gj V 

dx 7 


EEw«4» 

P,6 j,k 


where x = r/r a (r Q is a scaling radius). The parameter 


77 = 


Fq<? t 

4~r 0 c 


( 9 ) 


( 10 ) 


scales the effective optical depth of the flow (and is also related to the compactness, 
usually defined in terms of the luminosity, X, as L/r 0 = F 0 E/r 0 ). 

The experiments run for the non-expanding case, section 2 above, were repeated 
for the expanding case. Equations (9) for the spectrum were solved simultaneously 
with (5) and (7) for the flow. It was found that the geometry of the flow fairly quickly 
approached the usual adiabatic, mass conserving solution to (5) - namely, 7 ( 1 ") oc 
7 (r 0 )r/r 0 , for systems with only a moderate amount of pair creation/annihilation. It 
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was also found that emergent (frozen out) spectrum is predicted, given the creation 
spectrum, by the optical depth 

T(x ° )=r, L, iw (u) 

This is the quantity which measures the ’’number of collision times”, and can be directly 
compared to r = t/t 0 = cannot in the non-expanding runs. The emergent spectrum 
in a given case, that is for particular values of rj, -y 0 = 7 (zo) and a, particular creation 
spectrum, can be determined from the static solutions by choosing the scaled parameter, 
r, equal to r(x 0 ). For flows which are close to the adiabatic solution, ( 11 ) shows that 
t(x 0 ) — rj/ 7$ (since x 0 = 1 in general), so that the flow speed as well as the compactness 
determine the effective optical depth of the wind. 

5. Impact on the Surroundings and How the Flow is Observed 

The spectral evolution calculations were carried out in the co-moving (fluid) 
frame. The particle and photon spectra seen by an external observer will be Doppler 
shifted to higher energies. If the flow velocity at freeze-out, 7 i ast , is highly relativistic, 
then each quantum will be shifted up in energy by E ~ 7 lastE^ , if E? is that quantum’s 
energy in the fluid frame. 

The blueshifted photons, at least those with energies It several keV, axe likely 
to be observed directly. This calculation is thus roughly consistent with the observed 
hard X-ray/7-ray peak seen around hv ~ m e c ^ seen in AGX. 

The energetic pairs will not be seen directly. They can impact and interact with 
the local interstellar medium. For instance, they can act as the transrelativistic seed 
particles necessary for shock or turbulent Fermi acceleration to highly relativistic ener- 
gies. The particles can be seen indirectly through their radiation. Escaping positrons 
will of course annihilate on the local ISM; their high energies will result in a broad 
emissivity, centered ~ 1 - several MeV probably. Both species can produce synchrotron 
radiation, which would be expected to appear at radio to optical frequencies. 
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One might ask how synchrotron radiation can be apended to a calculation which 
was done for an unmagnetized plasma. The effect of a magnetic field would be to add 
synchrotron losses to the basic processes determining the spectral evolution, and to 
affect the flow dynamics. It was found that the first effect is unimportant for fields 
which satisfy 


B l « 


rrigC^ L cr e 


mfc 6 <7 C 77 


e 4 4zr 2 e e ^7 e 4 cr? r 0 e e 


( 12 ) 


if cr e is the cross section of the dominant reaction for the pairs, e e is the normalized 
mean energy of the pairs and equation (10) has been used. The magnetic fields are not 
important in the flow dynamics if 


B 2 «^8 : 

r z c ax r 


(13) 


where e is the mean (pair and photon) energy. For ” conventional” bright quasar num- 
bers, both (12) and (13) say that magnetic fields must be stronger than several hundred 
Gauss to have a strong effect on the evolution of the fireball. Fields weaker than this 
can still lead to detectable synchrotron emission without vitiating the unmagnetized 
calculation. 

Another possibly interesting effect of these escaping electrons was explored by 
the P.I. in collaboration with D. Baker (LANL/NASA), J. Borovsky (LANL) and G. 
Benford (U.C. Irvine). In this work, which was separate from the numerical calcula- 
tion above, it was noted that relativistic electrons might well be beamed, rather than 
isotropic, when they first encounter the nuclear (or jet) plasma in AGN. If this is 
the case, rapidly-growing electrostatic turbulence (such as Langmuir turbulence, from 
a streaming instability) can give rise to strong collective emission of electromagnetic 
waves. This phenomenon is observed in the lab, where it is found that this emis- 
sion can mimic the nonthermal, polarized emission usually interpreted as incoherent 
synchrotron. It was found that longer-wavelength electromagnetic instabilities could 
distort an initially ordered magnetic field and twist the collective emission into the ob- 
servers line of sight. Thus, it was argued that some part of the nuclear radio emission 
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observed from AGN might be due to collective processes rather than (or in addition 
to) incoherent synchrotron emission. 

Finally, an invited review was written on the role of high-energy particles and 
photons in the environment of massive black holes (which are thought to be the en- 
ergy source for AGN). Observations of high-energy quanta were reviewed, and general 
theoretical considerations as to the physics of these energetic quanta (both their initial 
’’creation” in nuclear accretion flows, and their ’’post-processing” according to the work 
described in this Final Report) were presented. 
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